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Abstract 

We investigate the correlation functions of the one-dimensional Asymmetric Simple 
Exclusion Process (ASEP) with open boundaries. The conditions for the boundaries 
are made most general. The correlation function is expressed in a multifold integral 
whose behavior we study in detail. We present a phase diagram of the correlation 
length. For the case the correlation length diverges, we further give the leading terms 
of the finite-size correction. 

1 Introduction 

As we all appreciate, exactly solvable models have been playing important roles for under- 
standing nonlinear physics in various areas. They make it possible to calculate physical 
quantities exactly so that we can enjoy arguments beyond approximation. What is more, 
the integrability of models often serves a deeper mathematical structure. For one thing, we 
consider the one-dimensional Asymmetric Simple Exclusion Process (ASEP). The ASEP 
is an exactly solvable model of stochastic process regarded as a primitive model of kinetics 
of biopolymerization, traffic flow, and so on. The model is known to have nonequilibrium 
properties in the stationary situation and is studied extensively for the nonequilibrium 
statistical mechanics |14| I15| I17| 1211 1221 I23| 120] . In the mathematical aspect, it has a 
connection with the theory of the (/-orthogonal polynomial |19[ H| l2*l|. 

The nonequilibrium statistical mechanics has attracted a renewed interest among physi- 
cists these days. Numerical experiments are more available than before and a lot of in- 
teresting applications are found out. Despite of these developments, the central principle 
of the nonequilibrium statistical mechanics is still unclear. In order to unveil the central 
principle, we need to look for basic concepts and clarify mathematical structures. There 
are some attempts ongoing for this end such as a construction of thermodynamics for the 
steady state which is far from equilibrium ^5j. At the same time, we may have recourse to 
reliable exactly solvable models and analyse the behaviors in detail. On the way, we can 
define quantities which are essential for the general theory. For example, for the ASEP we 
find a definition of the partition function which really is a generating function of the bulk 
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quantities as is familiar in the standard equilibrium statistical mechanics. In this sense, 
exact analytical arguments for individual models are of great value also in the study of 
the nonequilibrium statistical mechanics. 

Physically and mathematically, boundary conditions are important. In the case with 
open boundaries, an interesting feature of the ASEP is the connection with the Askey- 
Wilson polynomial. The Askey- Wilson polynomial is a one-variable ^-orthogonal poly- 
nomial. There is a set of lists called the Askey-scheme where many other orthogonal 
polynomials are recovered by parameter reduction from the Askey- Wilson polynomial 
PJIHIE]- As to physics, the Askey- Wilson polynomial also appears in various integrable 
models; the problem of Bloch electrons in a magnetic field on a lattice called the Azbel- 
Hofstadter problem which is closely related with the quantum Hall effect |25| . the XXZ 
model with open boundaries ^0]) the sine-Gordon model with a boundary and the more 
general algebraic formulation manipulating them (21 El- Through the connection with 
such a mathematical structure, we can obtain more information about the systems than 
we expect. 

The physics of the ASEP is rich. Because of the nonlinearity, the formation of shock 
profile happens [Sj. In fact, macroscopically the ASEP obeys the noisy Burgers equation. 
In the setting of an infinite system, the current fluctuation can be calculated based on 
the Random Matrix Theory through the statistics of the Young diagram [111 I16j . The 
Random Matrix Theory is a theory of universality classes. It is found that the current 
fluctuation of the ASEP belongs to the KPZ nonequilibrium universality class |12j . For a 
system with open boundaries, by the boundary effect there occurs phase transitions even 
in the one-dimensional system. In contrast, according to the Mermin- Wagner theorem, 
a phase transition does not take place in a one-dimensional equilibrium system. We can 
give a lot of fascinating phenomena like this. 

The aim of this paper is to study the correlation functions for the steady state of the 
ASEP with open boundaries. The ingredient of our analysis is the integral expression of 
the n-point function. The organization of the paper is as follows. In section 2, the model 
and related quantities of the ASEP are introduced. In section 3, based on j2lj, some known 
results are briefly summarized. It is remarked that the n-point function is expressed as 
a multifold integral. In section 4, the correlation functions for the steady state of the 
ASEP are studied and a phase diagram for the correlation length is presented. For the 
case where the correlation length diverges, we further inquire the integral representing the 
correlation function and obtain the leading terms of the finite-size correction. Section 5 is 
given for conclusion. 



Correlation function of ASEP with open boundaries 



3 



2 Asymmetric Simple Exclusion Process (ASEP) 

We introduce the model of the one-dimensional Asymmetric Simple Exclusion Process 
(ASEP). Roughly speaking, the model describes randomly hopping particles with exclusion 
interaction on a one-dimensional lattice. The term "asymmetric" refers to the asymmetric 
hopping of particles and "simple exclusion" refers to their hard-core exclusion interaction. 
We consider a finite system of size L with open boundaries at both left and right sides. 
On each site i = 1, • • • ,L of the lattice there can be only up to one particle, that is, if 
we denote by r< the particle number at site i, each site is either occupied by one particle 
Tj = 1 or empty n = 0. A particle hops to the nearest right (resp. left) site with the 
probability pndt (resp. pidt) during an infinitesimal time dt if the target site is empty. 
By rescaling time, we can set pr = 1 and pl = q- At the boundaries, a particle enters 
the system at the left (resp. right) boundary with rate a (resp. 5) if the site 1 (resp. L) 
is empty, and leaves at the left (resp. right) boundary with rate 7 (resp. j3) if the site 1 
(resp. L) is occupied. The boundary conditions are made most general containing four 
independent parameters a, (3, 7 and 5. Figure ^ illustrates the setup of the system. It is 
known that the system behaves very differently for the symmetric hopping case q = 1 and 
the asymmetric case q 7^ 1. We only consider the asymmetric case in the following. We 
can restrict q < 1 because particle-hole exchange simply gives q > 1. 
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Figure 1. The ASEP with open boundaries. Black circles and boxes represent particles and 
lattice sites, respectively. 

The state of this system is labelled by the configuration of particles {t^} = (n, T2, ■ • • , tl) 
where t% = 1,0. Then, the time evolution of the probability distribution P({ri}; t) is gov- 
erned by the master equation 



d_ 
dt 



P{{r i };t)) = H\P({T i };t)). 



(2.1) 



Here we used the bra-ket notation. The time evolution operator H is a matrix whose 
elements are the transition probability rates, and is defined explicitly by 



^V+i + oT<a " ~(1 " of )(1 + aUt) - |(1 + of )(1 - af +1 ) 



+ aat + 7 ^r " 2 + a ^ ~ ^ ~ °l) + 6(7 1 + P°L ~ ^ + °i)-%0-- O 

(2.2) 

where erf (a = x, y, z) is the Pauli matrix acting on the i-th site in the basis T{ = 1,0 and 
af = (erf ± \J — laf )/2. From the general properties of the time evolution operator, there 
exists a unique steady state: the state which is stationary under the time evolution, or the 
highest eigenstate of H with eigenvalue zero. In this paper, we focus on the steady state. 
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As is obvious from the explicit form of the time evolution operator, the one-dimensional 
ASEP can be regarded as a non-hermitian generalization of spin chain models in which the 
spin-flip occurs asymmetrically. There are a lot of techniques to solve spin chain models. 
For example the Bethe ansatz method is the representative one. For our purpose, the 
matrix method is useful [ZJ. According to the matrix method, the probability distribution 
of the steady state which we denote by P({rj}) is expressed in the matrix product form: 

L 

P{{n}) = ±{W\ + (1 - n)E)\V). (2.3) 

ZL f=l 

Here D and E are matrices; D (resp. E) corresponds to the occupation (resp. emptiness) 
of a particle. (W\ and |V) are vectors; (W\ (resp. \V)) corresponds to the left (resp. 
right) boundary. One can show that Q2.3|) gives the exact steady state of the system if the 
matrices and the vectors satisfy [7] 

DE - qED = D + E, (2.4) 
(W\(aE - jD) = (W\, (2.5) 
(PD-5E)\V) = \V). (2.6) 

The normalization constant Zl in (|2.3|) is given by 

Z L = (W\C L \V) (2.7) 

where C = D + E. We call Zl the partition function and C the transfer matrix. The 
reasons are as follows. In the similar way to the standard statistical mechanics, the bulk 
quantities of the system can be calculated from Zl- More precisely, if we introduce the 
fugacity A and define Zl(X) = (W\(XD + E) L \V), the density p and the density fluctuation 
Ap 2 are obtained from 



p = ±. x JL logZLix) 



(2.S 

A=l 



Ap 2 = log Z L (\) 



(2.9) 

A=l 



The definition of the matrix C implies that it has all information on the possibility at each 
site, occupation and emptiness of a particle. Since it is also site-independent, the partition 
function is calculated from the eigenvectors of C. Furthermore, the partition function is 
dominated by the contribution of the largest eigenvalue of C in the thermodynamic limit 
L — > oo. Interestingly, phase transitions occur when the set of eigenvalues of C varies 
qualitatively. We shall come back to this point later. 

Among other physical quantities we are interested in the n-point functions. For exam- 
ple, the one-point function and the two-point function are respectively written as 

{n) = ^-{W\C- 1 DC L - i \V), (2.10) 
Zl 

( TiTj ) = ^-{Wp^DG-^DC^^V). (2.11) 
Zl 
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Here the bracket in the LHS means the average over the steady state distribution (j2.3j) . 
An important nonequilibrium quantity, the particle current J, is defined by the two-point 
function: J = (rj(l — Tj+i) — q(l — T^Tj+i). In terms of the matrix product expression, we 
have 

J=^i. (2.12) 

Indeed, the current J takes nonzero value even in the steady state. Therefore, the ASEP 
is considered to be a nonequilibrium model. 



3 Known results 

We briefly review known results for the steady state of the ASEP with open boundaries. 
Most of them are found in [21]. There, the Askey- Wilson polynomial plays a chief role. The 
Askey- Wilson polynomial P n (x) = P n (x; a,b, c, d\q) ,n £ Z + is a g-orthogonal polynomial 
with four parameters a, b, c, d besides q. It is regarded as the most general one in the 
hierarchy of the one- variable (/-orthogonal polynomial family in the Askey scheme [T|l?Hll3|. 
The orthogonality relation for P n (x) is 



I I p m ( x )p n ( x )dx + V V w f k P m (x{)P n (x f k ) = h n 5 r , 



°n' J rnn 

f=a,b,c,dk:K\fq k \<\f\ 

(3.1) 

where, by using the notation of the g-shifted factorial (oi, • • • , a s ; q) n = Y\l=i IIfc=o(l ~~ 

a r q k ), 



e 2ie ,e- 2id ;q) 



h n (1 — q r ' abcd){q, ab, ac, ad, be, bd, cd; q) n 
ho (1 — q 211 " 1 abed) (abed; q) n 

(abed; c 
, ab, ac, ad, be 



(3.3) 



ho = (abed; q)^ ^ 

(q, ab, ac, ad, be, bd, cd; q)oo ' 

(3.5) 



a (a ,g)oo (1 - a 2 q 2k )(a 2 ,ab, ac, ad; q) k / q y 

k (q, ab, ac, ad, a~ 1 b, a~ l c, a~ l d; q)oo (1 — a 2 )(q, ab~ l q, ac~ 1 q, ad~ l q; q)oo \abcd) 

(3.6) 

and similarly for f = b,c, d. In the theory of the orthogonal polynomial, the three-term 
recurrence relation and the orthogonality relation are fundamental objects. The three- 
term recurrence relation is represented by an infinite-dimensional matrix called the Jacobi 
operator T as an eigenvalue equation, 

T\P{x)) = x\P{x)} (3.7) 

where \P(x)} = t (Po(x), Pi(x), ■ ■ ■)■ The fact is known as the spectral theorem that the 
spectrum of the Jacobi operator corresponds to the support of the orthogonal measure of 
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the orthogonal polynomial (see for example [H]). From (|3.1j) we can read off the spectrum of 



T as [—1,1] U{x[,f = a,b,c,d\l < \fq k \ < \f\,k G Z + } for the Askey-Wilson polynomial. 

For the steady state of the ASEP, we found in [21] that the transfer matrix C is directly 
related to the Jacobi operator of the Askey-Wilson polynomial T: 



C = 2 + 2T 
with parameters being set 

a = K a n : b 

where 



k Qj7 , d — 



a,7 



1 

2a 
1 



(1 



a 



+ 7) ± a/(1 _ 1 - a + 7) 2 + 4cry 



(1 - g - /? + 5) ± vXL - q - 13 + 5) 2 + 4/3<5 



(3.8) 
(3.9) 

(3.10) 
(3.11) 



Therefore, the eigenvector \x) of C is written in terms of the Askey-Wilson polynomials 
and the spectrum of C is determined as [0, 4] U {2 + 2x^ f = a, fell < \ fq k \ < |/|, k € Z + }, 
noticing that |c|, |d| < 1. In fact the partition function (J2.7|) is calculated to be in the 
integral form, 



dz (^z-2. g)oo[(1 + z)(1 + z -l )/(1 _ g)] i 

q 'iiriz (az, a/z, bz, b/z, cz, c/z, dz, d/z; q)^ 



(3.12) 



where the integral contour C is a closed path which encloses the poles at z = aq k , bq k , cq k , 
dq k (k £ Z + ) and excludes the poles at z = (og fe ) _1 , {bq k )~ l , (cq k )~ 1 , (dq k y 1 (k G Z + ). 
This integral is actually the moment integral with respect to the weight function of the 
Askey-Wilson polynomial. 

Although the integral of (|3.12|) is difficult to evaluate, physically important thermody- 
namic limit L — > 00 makes it possible to obtain the leading contribution. In this limit, the 
largest eigenvalue of C dominates the contribution to the partition function. Indeed, it is 
clear if we expand 



(W\C L \V) =Y,(W\C L \x)(x\V) 

X 

= ^2(W\x)x L (x\V) 



— const. x msx . 

Taking account of the largest eigenvalue, we find that there are three phases. The bulk 
quantities for each phase are eventually obtained as the following. 
(A) a > 1 and a > b: x max = (1 + a)(l + a -1 ), 



Z L ~ 
J ~ (1 



(q,ab,ac,ad,a 1 b,a 1 c, a 1 d;> 

a 1 

p ~ 



(l + o)= 



1 + a 



(1 + aXl + a- 1 ; 



Ap 2 ~ 



(l + a) 2 L' 



(3.13) 
(3.14) 
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(B) b > 1 and b > a: x max = (! + &)(! + b~ x ), 



Z L 



(b- 2 ;q) c 



(1 + &X1 + 6- 1 ) 



(q,ba,bc,bd,b 1 a,b 1 c,b 1 d;q) OQ [_ 1 



(1 + 6)2 



1 + 6' 



(1 + 6) 2 L' 



(3.15) 
(3.16) 



(C) a < 1 and 6 < 1: x max = 4, 



J 



(^; q) 



1, c, d; ^ \/vri 3/2 V 1 - Q 



l-q 



(3.17) 
(3.18) 



Note that a phase transition occurs when the largest eigenvalue x max of the transfer matrix 
C changes from one to the other. 

Now, let us consider the correlation functions. We have the following integral formula 
for the n-point function, 



'Jn,/ 



lim 



(1 -q) L Z L fr,- 



dz r 



n+l 

n 7 \-i: n , 



-1 



71+1 

H (4,z" 2 ;( ? ) 00 [(l + z m )(l + l/^)]^-^- 1 

m=l 

(azi , a/21, czi , c/zi , 6z n+ i , 6/ z n+ i , dz n+ i , d/ ; q) c 
1 



n 

m=l 



"^(^-m+1 ^J7i ) 



+ 



(q; g) 2 x3 (z m+ i + l/z m+ i - j m z m - j m j ' z m ) 

(^m^m^m+l i Cm^m/ Zm+li / '^mi ^m/^m^m+li (?)c 



(3.19) 



where jo = 0,j n+ \ = L + 1 in the product. Regularization parameters £ m 's are included 
since otherwise this integral becomes a singular integral. The contour C m (m = 2, • • • , n + 
1) encloses the poles at z m = £,m~iz m -iq k , ^m-iz'^_ 1 q k (k G Z + ) and excludes the poles 
at z m = ((,m-iZm-iq k )~ 1 j (^m-i-^m 1 -!^) 1 G ^+)- I n addition, for m = 1, C\ encloses 
the poles at zi = aq k ,cq k (k G Z + ) and excludes the poles at z\ = (ag fc ) _1 , (cg fc ) _1 
(k G Z + ). Also, C n+ i encloses the poles at z n+ \ = bq k ,dq k (k G Z+) and excludes the 
poles at z n +i = (bq k )~ 1 , (dq k ) _1 (k G Z + ). In the thermodynamic limit, the second 
largest eigenvalue also comes into the game. For the two-point function (|2,11|) . we have 
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an expansion 

— "^max (^^l-^max) (^max|-C|2?max) (^maxl^l^max) (^maxl^") 



max | D | X max ) (x max | V) 

\ 3?max , 



( ^2nd 



max 



( X2n& 



x max (W\x max ){x max \D\x 2 nd){x2nd\D\x mlLX ){x max \V) 



L-j 



^max max|-^|^max) (^max \D\x 2n d)(x2nd\V). 

\ ^max / 

From this expansion, it is found that the ratio a^nd/^max determines the decay scale of the 
correlation. However, the second largest eigenvalue is defined only for the phases (A) and 
(B) since the spectrum of C is continuous for the phase (C). For this reason we investigate 
in detail the correlation functions in the thermodynamic limit in the subsequent section. 



4 Correlation functions 

We will study the n-point function. In the following, suppose that the site numbers are of 
order or suborder L, and take a thermodynamic limit L — > oo. 

First, the phases (A) and (B). In these cases, the second largest eigenvalue of C is 
separated from the largest one and can be explicitly read off from the spectrum. It 
is not necessary to touch the integral (|3.19|) directly. The correlation functions decay 
exponentially and the correlation length is determined by £ = (hiXmax/a^nd) -1 - The 
phase (A) is further classified into three sub-phases. 
(Ai) aq > 1, aq > b, 

Xmax = (l + aXl + cr 1 ), x 2n d = {l + aq)(l + (aq)- 1 ). (4.1) 
(A 2 ) a > b> aq, b > 1, 

Xmax = (1 + aXl + a- 1 ), x 2nd = (l + 6)(l + 6- 1 ). (4.2) 
(A3) a > 1 > aq, b < 1, 

Xmax = (l + a)(l+a _1 ), x 2nd = 4. (4.3) 

Because of the particle- hole symmetry, the phase (B) is classified in the same way except 
a <-> b in the above. 

Next, the phase (C). As mentioned before, since the spectrum of the transfer matrix C 
is continuous, the second largest eigenvalue is not defined. Thus, we should take care of the 
integral representation. The behavior of the correlation function turns out to be something 
different from the exponential decay. We can identify the correlation length with infinity 
in this case, and finally we obtain a phase diagram with respect to the correlation length 
as Figure [5J The horizontal (resp. vertical) axis corresponds to the left (resp. right) 
boundary. As 1/a increases, the input from the left boundary tends to increase; as 1/6 



Correlation function of ASEP with open boundaries 



9 




1/a 



Figure 2. Phase diagram for the correlation length. The solid lines distinguish the three phases 
(A), (B) and (C). The dashed lines separate the three sub-phases in (A) and (B). 



increases, the output into the right boundary tends to increase. The special case 7 = 5 = 
was argued and a similar phase diagram was obtained in |19j . 

We consider the finite-size correction of the correlation functions. Taking a close look 
at the integral IjM.mj) , we can verify a decomposition 



k=0 {h,- ,i k }c{ji,- ,j n } 



(4.4) 



By definition, (St^ ■ ■ ■ Sri k ) is the term represented by some integral which depends only 
on {ii, ■ ■ ■ , ifc} out of {ji, ■ ■ ■ ,j n } and is independent of n. These terms come from the 
existence of the boundaries. Therefore, we can identify them as the finite-size corrections. 
The n-point correction term of the correlation functions is written in the following integral 
form, 



( St ji • • • Stj J 



n+l 

n 

,m=l 
n 



2n 

oo 



(1 - q) L Z L 



2vr 



n+l 

IK 

m=l 



2id m e — 2 i6 m \ 

! c iH)oo 



(2 + 2 cos e r , 



\jm jm — 1 1 



n 



{2cos6 jm+1 - 2cos0 im ) 



- ( i(6i m +9i ,,) i(9i„,-9j 



(4.5) 



where jo = 0, j n +i = L + 1 in the product. In the thermodynamic limit, the contribution 
to this integral is concentrated on 9 m = 0,m = ,n + 1 due to the presence of 
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(2+2 cos 9fn)^ m "'™ 1 • Approximating (|4.5j) around this point by the saddle point method 
and substituting the result Q3.17[) . we arrive at a simpler expression, 



n+l 



7T 



(2tt; 



n+l 

n 

m=l 



II #ni ex P 



d0, 



m=l 



1 



|Q (#m - #m+l) 



m=l 



(4.6) 



The integrands have singularities at # m = # m +i, but this multifold integral should have a 
certain finite value because we are evaluating a finite-valued physical quantity. To clarify 
this point, we rewrite the integral in the RHS; in fact it is equivalent to another form, 

2n(£l, ' ' ' > £n+l) 

n+l "+1 



n 



o 



Em £ ^^m — 1 



1 



Z^fc=l sfc Z^fc=l '/A 
(4.7) 

with £ m = (j m — jm-i — l)/4, m = 1, • • • , n + 1. In turn, the integral in (|4.7|) is integrable 
and definitely has a finite value. Rewriting of the integral is performed as follows. First 
of all, we have the Gaussian integral 



2y^ 



Then, integrating the above over r\ from to £ yields 



°° e-W 2 -* 2 ) 
9 2 -ip 2 



-d9 + 



dO 



9 2 



2y^ 



The second term of the LHS is equal to ^-S(tp). Now multiplying e ^ on both sides and 
then differentiating by £, we have 



o # 2 -^ 



-& 2 



lo 2^? 

Note that we can drop ip 2 e~^ lp2 S(ip). Applying this formula to the integral of 9\, then #2, 
■ ■ ■ , n , we get 

n+l 

II #m exp [-infil 



n+l 

n 

m=l 



? 2 1 



d0 r 



m=l 



n (oi - 9 



2 i 

m+li 



m=l 



a 



n £)n 

n=l x 



6» n+1 exp 



m=l 
/n+l 



Em £ sr^rn — 1 



vfc=l 



fc=l 



d0 



n+l- 
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Finally, performing the 9 n +i integral, we obtain (|4.7|) . After all, we have the following 
result. 



r 3/2 V 71- t I ii - 1 h 
(2^)™+i " 



j n -l - 1 L-j r 



(4.8) 



From this expression, it is found that the correlation functions are independent of the 
boundary parameters a, b, c, d and q. Moreover, one can show that the n-point correction 
(Srj 1 ■ ■ ■ Srj n ) is of order L~ n l 2 if ji, ■ ■ ■ ,j n are of order L. 

We obtain exact expressions for the n-point corrections by calculating (|4.7|) explicitly. 
We show the first few examples: 

1 i-L/2 



(5Ti) ~ 
{STiSTj} 



lM 2 ^^/i{L-if 
L 



(4.9) 



1 

27T 



-i/2 VT^iH 2 + U ~ i? + ijL - j) — (L — jf) 

2y/i(L-j)j(L-i) 



r -i3 i( L ~j) 
+L - arctan \ / ; — 

2 V U - { ) L 



(4.10) 



(6Ti5Tj8T k 



1 



2W 2 



L -l/2 



2^/i(L-k)j(L-j)k(L-i) 



x{i 3 + i 2 (j -i)- i(j - if - (j - if - i(j -i)(k- j) - (j - if(k - j) 
+U ~i)(k- jf + (k — jf -i 2 (L-k)- i(j -i)(L-k) + i(k - j)(L - k) 



+{j-i){k-j){L-k) + {k- 
(L — kf + 6k(L — k) 



+L -3/2 
_ L -3/2 



jf(L -k)+ i(L - kf -(k- j)(L -k)-{L- kf} 
- 3k 2 



2\/Z~ 



i 2 + 6i(L 



kk 3 / 2 

- 3(L - if 



■ arctan 



l i(k-j) 

U -i)k 



arctan , 



(j-i)(L-k) 



(4.11) 



2V~i{L-if/ 2 \{k-j){L-i) 

Physically, fluctuation of the correlation function is related to the connected part, 

• • • ^iJconn = < {5r h - (5T n )) ■ ■ ■ {St 5i - (<%)) >. (4.12) 

Figure|3]iHustrates the behavior of the connected three-point correction term (dTidTjSrk) conn 
for large L. It is interesting that there can be a maximal or minimal point for the par- 
ticle occupation according to the positions of the other particles. One also notices that 
(dTidTjdrk) conn vanishes at the boundary. This is a physical consequence of the fact that 
particles are in equilibrium with the reservoir at the boundary. We can show in general 
that the connected n-point functions vanish at the boundaries. From (J4.7j) . we directly 
have 



1. 



1 f^ 1 dt] d 

V£i Jo Vv <9£i 
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Figure 3. The connected part (STiSTjSrk) CO im for large L. The two curves are plotted as a 
function of the position of the leftmost particle i with the position of the middle particle j and 
that of the rightmost one k fixed at the place of the small black circles and the big black circle 
respectively. The solid line is for j = 0.1L and k = 0.5L. The dashed line is for j = 0AL and 
k = 0.5L. 



Note if we take £i —* 0, only the first term of RHS is singular and the second term vanishes. 
Then, we have lii)ij_»o (<$Ti#Tj • • • ^ T fc) — (^ T i)(^ T j ' ' ' ^ T k) = because (Sri) ~ ^7^7 i — > 0, 

and furthermore limj_>o 0>Ti(8Tj — {Srj)) ■ ■ ■ {St^ — {St^))) — (<5ri)^(<5rj — {6tj)) ■ ■ ■ (5rk — 
(<5Tfc))y = , which establishes 

lim {5T h ■ ■ ■ 5T jn ) conn = 0. (4.13) 

Similarly, (5tj 1 ■ • • 5Tj n ) conn vanishes also in the limit j n — > L. 

What we have calculated here are the leading terms of the finite-size corrections of 
(5rj 1 ■ ■ ■ 5Tj n ) conn . It is also possible to have 1/L-expansion of it by elaborating approxi- 
mations of diSJ and (|3~T7|) . 



5 Conclusion 

We have investigated the correlation functions for the steady state of the ASEP with open 
boundaries in the thermodynamic limit. Main results are summarized as follows: 1) a 
phase diagram is given with respect to the correlation length £, 2) for the case £ = oo, the 
leading terms of the finite-size corrections are calculated accurately. 

For the future problem, we are interested in the time-correlations of the ASEP. In order 
to calculate the time-correlations, one should know all the eigenvalues and eigenvectors 
of the time evolution operator H. At present, we have an exact information only for the 
steady state. It would be possible to solve the eigenvalue problem for H by means of the 
Bethe ansatz method and evaluate the time-correlation functions for the ASEP with open 
boundaries. 
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In the nonequilibrium statistical mechanics, we like to have the general theory for the 
relations among nonequilibrium quantities. For the purpose, the study of the current 
fluctuation is primarily important. An attempt to calculate the current fluctuations of 
the ASEP with open boundaries was made for the symmetric case q = 1 6 . There, the 
authors considered a hierarchy of equations for the correlations without directly solving 
the eigenvalue problem for H. Whether their approach is applicable to the asymmetric 
case q ^ 1 is one of open questions. 
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